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The charged-current quasi-elastic scattering of muon neutrino on the oxygen target is 



in , 

•^ ' analyzed for neutrino energy up to 2.5 GeV using the Relativistic Distorted- Wave Impulse 

•/^ ' Approximation (RDWIA). The inclusive cross sections (P'a/dQ^ ^ calculated within the RD- 

WIA, are lower than the Relativistic Fermi Gas Model (RFGM) results in the range of the 
square of four-momentum transfer Q^ <0.2 (GeV/c)^. We have also studied the nuclear- 
^ model dependence of the neutrino energy reconstruction accuracy using the charged-current 

^ ' quasi-elastic events with no detector effects and background. We found that for one-track 



events the accuracy is nuclear-model dependent for neutrino energy up to 2.5 GeV. 



^ . PACS numbers: 25.30.-c, 25.30.Pt, 13.15.4-g 

(N 
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I. INTRODUCTION 

o 

QQ ' The field of neutrino oscillations is rapidly developing. The goals of the current and planed 



set of accelerator-based neutrino experiments [1-7] are the precision measurements of observed 
neutrino mass splitting and the detailed study of the neutrino mixing matrix. The data of these 



j^ ■ experiments will greatly extend the statistics due to extremely intense neutrino beamline. 

To study the neutrino oscillation effects on the terrestrial distance scale, the neutrino beams 
cover the energy range from a few hundred MeV to several GeV. In this energy range, the dominant 
contribution to the neutrino-nucleus cross section comes from the charged-current (CC) quasi- 
elastic (QE) reactions and resonance production processes. The cross section data in this energy 
range are rather scarce and were taken on the targets, which are not used in the neutrino oscillation 
experiments (i.e. water, iron, lead or plastic). 

In this situation, the statistical uncertainties should be negligible as compared to systematic 
errors in the incident neutrino flux, neutrino interaction model and the detector effects on the neu- 
trino events selection and neutrino energy reconstruction. Apparently, these uncertainties produce 
systematic errors in the extraction of oscillation parameters. 

To evaluate the neutrino mass squared difference in the muon neutrino disappearance exper- 



iments, the probability of f^ disappearance versus neutrino energy is measured. Because the 
CCQE interaction represents a two-particle scattering process, it forms a good signal sample, and 
the neutrino energy may be estimated using the kinematics of this reaction. There are two ways 
of measuring the neutrino energy: either kinematic or calorimetric reconstruction. In detectors 
with the energy threshold for proton detection e^^ > 1 GeV (Cherenkov detectors) the muon neu- 
trino CCQE interactions will produce the one-track events, i.e. only muons are detected in the final 
states. The kinematic reconstruction is applied for these events. Assuming the target nucleon to be 
at rest inside the nucleus, the correlation between the incident neutrino energy and a reconstructed 
muon momentum and scattering angle is used in this method. 

In the fine-grained detectors the two-track CCQE events are detected, and the calorimetric 
reconstruction can be applied, if the particle identification of the second track and the resolution 
for proton momentum are reliable. In this case the visible neutrino energy is simply a sum of the 
reconstructed muon energy and kinematic proton energy. In this paper we consider the procedures 
for neutrino energy reconstruction, which are based on the kinematics of the CCQE interaction. 

In general, the detector efficiency and energy response are highly dependent on the type of 
interaction: QE or non-QE (the resonance and deep inelastic scattering). The Monte Carlo (MC) 
simulation of the detector response to neutrino interactions is used for tuning the values of cuts 
for separation of the QE and non-QE (nQE) events and for estimating the efficiency of detecting 
the CCQE events after all cuts. To model the scattering from a nuclei, the most part of an event 
generator [8] is based on the RFGM [9], in which the nucleus is described as a system of quasi-free 
nuclei with a flat nucleon momentum distribution up to the same Fermi momentum pp and nuclear 
binding energy ei,. But this model does not take into account the nuclear shall structure, the final 
state interaction (FSI) between the outgoing nucleon and residual nucleus and the presence of 
short-range nucleon-nucleon (NN) correlations, leading to appearance of a high-momentum and 
high-energy component in the nucleon momentum-energy distribution in the target. 

The comparison with the high-precision electron scattering data has shown [10] that the accu- 
racy of the RFGM prediction becomes poor at low Q^, where the nuclear effects are largest, and 
this model fails [13] in application to exclusive cross sections. The modern quasi-elastic neutrino 
scattering data (the CCQE event distribution as a function of Q^) [2, 11, 12] also reveal the inad- 
equacies in the present neutrino cross section simulation. The data/MC disagreement shows the 
data deficit in the low Q^ region (Q^ < 0.2 GeV^) and the data excess in the high Q^ region. The 
disagreement at low Q^ would eventually result in the data/MC disagreement in the reconstructed 
neutrino energy. 



The Relativistic Distorted- Wave Impulse Approximation, which takes into account the nuclear 
shall structure and FSI effects, was developed for description of electron-nucleus scattering, and it 
was successfully tested against the data [14]. The RDWIA approach was also applied to neutrino- 
nucleus (lyA) interactions for calculating the exclusive and inclusive QE cross sections [13, 15-17]. 
In Ref. [13] the FSI effects on the inclusive cross section in the presence of the NN-correlations 
were estimated. 

The aim of this work is twofold. First, we compute the RDWIA CCQE cross section versus 
Q^ for muon neutrino scattering off oxygen. Second, we show the nuclear-model dependence of 
the efficiency of two-track CCQE events selection. We also estimate systematic uncertainties in 
the reconstructed neutrino energy within the RDWIA and RFGM taking into account the nucleon 
momentum distribution in the target, i.e. the nucleon Fermi motion effect. 

The outline of this paper is as follows. In Sec. II we present briefly the formalism for CCQE 
scattering process and the RDWIA model. The nuclear-model dependence of cuts, which are 
applied for CCQE events selection, as well as the neutrino energy reconstruction methods are 
discussed in Sec. III. The results of numerical calculations are presented in Sec. IV. Our conclusions 
are summarized in Sec. V. In Appendix A we present the equation for neutrino energy, and in 
Appendix B the expressions for the moments of the reconstructed neutrino energy distribution are 
given. 

II. FORMALISM OF QUASI-ELASTIC SCATTERING AND MODELS 

We consider the electron and neutrino charged-current QE exclusive 

u{ki) + A{pa) -^ l{kf) + N{p,) + B{pb), (1) 

and inclusive 

u{ki) + A{pa) -^ l{kf) + X (2) 

scattering-off nuclei in the one-W-boson exchange approximation. Here / represents the scattered 
lepton (electron or muon), ki = (cj, fcj) and kf = (sf, kf) are the initial and final lepton momenta, 
PA = {£a-,Pa)-> ^^*^ Pb = {^b,Pb) ^^6 the initial and final target momenta, Px = {ex,Px) is the 
ejectile nucleon momentum, q = (lo, q) is the momentum transfer carried by the virtual W-boson, 
and Q^ = —q^ = q^ — uP' is the W-boson virtuality. 



A. CCQE neutrino-nucleus cross sections 

In the laboratory frame the differential cross section for the exclusive (anti)neutrino CCQE 
reaction, in which only a single discrete state or narrow resonance of the target is excited, can be 
written as 

rfV -r\p^\M91^^11il,^w^^, (3) 



dsfdUfdn^ (27r)5 Si 2 

where Qf is the solid angle for the lepton momentum, 0,x is the solid angle for the ejectile nucleon 



momentum, G ~ 1.16639 xlO MeV is the Fermi constant, 9c is the Cabbibo angle (cos^. 



c 



0.9749), Lfj_i, and W^"^ are, respectively, the lepton and weak CC nuclear tensors. The recoil factor 
R is given by 

R= dexS^Ex +eB -(jJ - niA) 

and Ex is the solution to the equation 



-1 
1 _ ^x Px • Pb 



^B Px ■ Px 



(4) 



Ex+eb -ruA-uj = 0, (5) 

where eb = \/"^b +Pb' Pb = Q ~ Px^ Px — V^x ~ "^^; ^i^d niA, ms, and m are masses of the 
target, recoil nucleus and nucleon, respectively. The missing momentum pm and missing energy 
Em are defined by 

Pm=Px-<l (6a) 

Em =m + niB - ruA (6b) 

The lepton tensor can be written as a sum of symmetric L^ and antisymmetric L^ tensors 

L^''' = Lf + L^ (7a) 



Lf = 2 (k^k) + k\e^ - g'^^'kikfj (7b) 

L>'/ = h2te>^''^^{hUkf)p, (7c) 

where h is +1 for a positive lepton helicity, and —1 for a negative lepton helicity, e'^'^'^^ is the 

antisymmetric tensor with e0i23 _ _^^^^^ _ ^ The weak CC hadronic tensors W^i, are given by 

bilinear products of the transition matrix elements of the nuclear CC operator J^ between the 
initial nucleus state \A) and the final state \Bf) as 

W^, = Y,{Bf,Px\JM){A\4\Bf,Px), (8) 
/ 



where the sum is taken over undetected states. 

In the inclusive reactions (2) only the outgoing lepton is detected, and the differential cross 

sections can be written as 

dV 1 \kf\ G''zosH 

dsfd^f ~ (27r)2 e, 2 ''^"^^ ' ^^> 

where W^^ is an inclusive hadronic tensor. In the reference frame, in which the z axis is parallel 

to the momentum transfer q = ki — kf and the y axis is parallel to ki x kf, the exclusive neutrino 

scattering cross sections take the forms 



d^(^ \Px\^x^2 



G cos 6c£f\kf\R{^VQRQ + vtRt + vttRttcos2(J) + VzzRz 



d£fdQ.fd^x (2vr) 

+ (vxzRxz - voxRox) cos 4> - vqzRqz + h [vyz (Ry^ siu (j) + Ry^ cos (p) 

- VOyiR'oy Sm<j) + Roy cos (p) - V^yRxy] }, (10) 

where 6, ip are lepton scattering angles, and 9x, (p ^^^ outgoing nucleon angles, Vi are the neutrino 
coupling coefficients, and Ri are independent response functions [13], which depend on the variables 
Q^, to, \Px\, and O^- Similarly, the inclusive lepton scattering cross sections are reduced to 

, ,^ = — — 2-^e/|fc/|(foi?0 + VtRt + VzzRzz - VOzRoz " hv^yRxy) (H) 

where the response functions now depend on Q^ and uj only [13]. 

We describe the lepton-nucleon scattering in the Impulse Approximation (lA), in which only 
one nucleon of the target is involved in the reaction, and the nuclear current is written as a sum 
of single-nucleon currents. Then, the nuclear matrix element in Eq.(8) takes the form 

{p,B\J^\A) = f d^r exp{it-r)'^^'\p,r)Tt'<^{r), (12) 

where F^ is the vertex function, t = ebq/W is the recoil-corrected momentum transfer, W = 
\/{mA + ^)^ ~ Q^ is the invariant mass, <5 and '^^' are relativistic bound-state and outgoing 
wave functions. 

The single-nucleon charged current has V—A structure J'^ = Jy + Jj^. For the free- nucleon 
vertex function F'^ = Ty + F^;^ we use the CC2 vector current vertex function [18] 



F^ = Fv{Q^)l^ + i^'^'^FMiQ') (13) 



^1^ = F,.('n2W^ ^ irrl^'^JilLF.Jn^ 

and the axial current vertex function 

r^ = FAiQ^h^-/5 +P {Q^)q^l5, (14) 



where a^'^ = i[y^,'j'^]/2. The weak vector form factors Fy and Fm are related with corresponding 

electromagnetic 

current (CVC) 



(eO (eO 

electromagnetic factors for proton -Fj „ and neutron -Fj„ by the hypothesis of conserved vector 



p. — p(eO _ p(eO 
* i,p i,n ' 



(15) 



where Fy and F^ are the Dirac and Pauli nucleon form factors. We use the approximation 
of [19] for these form factors. Because the bound nucleons are off shell, we employ the de Forest 
prescription for off-shell vertex [18] and the Coulomb gauge - for vector current Jy- 

The axial Fj^ and psevdoscalar Fp form factors in the dipole approximation are parameterized 
as 



Fa{Q' 



Fa{0) 



FpiQ' 



2mFAiQ^ 



(1 + Q2/Ml)2' "^^ ' m2+Q2 ' 
where -Fa(O) = 1.267, Ma — 1.032 GeV is the axial mass, and m-,^ is the pion mass 



(16) 



B. Models 



In the independent particle shell model the relativistic bound-state functions $ in Eq.(12) are 
obtained within the Hartree-Bogolioubov approximation in the a — to model [20]. The bound-state 
spinor takes the form 



^Km(r) 






(17) 



where 



u nis m 



(18) 



is the spin spherical harmonic, and the orbital and total angular momenta are given, respectively, 
by 

1. 1 



SJk + 



T 2 



3 I^K^^ 



1 



(19a) 
(19b) 



with Sk = signi^K). The missing momentum distribution is then 



"m)\ I , 



(20) 



where 

^k{p) = dr r'^je{Pmr)F^{r) (21a) 

G.^ip) = I drr^je{Pmr)G.^{r), (21b) 



and j^{x) is the Bessel function of order ^, and H.' = 2j — L In this work the bound- nucleon wave 
functions [21] are used in the numerical analysis with the normalization factors Sa relative to 
full occupancy of ^^O: ^'(Ips/a) = 0.66, S'(lpi/2) = 0.7 [14] and ^(Isi/a) = 1. Note that the 
calculation of the bound-nucleon wave function for 1^3 /2-state includes the incoherent contribution 
of the unresolved 2si/2C^5/2 doublet [22]. 

In the RDWIA the ejectile wave function ^ in Eq.(12) is obtained following the direct Pauli 
reduction method [23, 24]. It is well known that the Dirac spinor 



* = (22) 




can be written in terms of its positive energy component ^_|- as 

\e+M+S-V^ + '/ 
where S = S{r) and V = V{r) are the scalar and vector potentials for the nucleon with energy E. 
The upper component ^-|- can be related to the Schrodinger-like wave function ^ by the Darwin 
factor D{r), i.e. 



<ll+ = ,/W)i. (24) 

^ E^M + S(r)-Vir) 

^ ' E + M ^ ' 

The two-component wave function ^ is the solution of the Schrodinger equation containing equiv- 
alent central and spin-orbit potentials, which are functions of the scalar and vector potentials S 
and V , and are energy dependent. We use the LEA program [25] for numerical calculation of the 
distorted wave functions with EDADl SV relativistic optical potential [26]. 

In the Plane- Wave Impulse Approximation (PWIA) the final state interaction between the 
outgoing nucleon and the residual nucleus is neglected, and the nonrelativistic PWIA exclusive 
cross section has a factorized form [27] 



dejdflfdfl^ 



Kae^V{E,p) (26) 



where K = i?p^e^/(27r)^ is the phase-space factor and a^x is the half-off-sheh cross section for 
neutrino scattering by a moving nucleon. The nuclear spectral function V{E,p) can be written as 

V{E,p) = Y,\ {Bf\a{p)\A) \^5{E-e^) (27) 

/ 

and the nucleon momentum distribution Pf^ip) for the orbit /? is related to the upper component 

of the corresponding bound-state wave function (21a) as 

2 

27r2 



P^ip) = p- F^ip) . (28) 



According to the JLab data [14], the occupancy of the independent particle shell model orbitals 
of -"^^O equals about 75%, on the average. In this work we assume that the missing strength (25%) 
can be attributed to the short-range NN-correlations in the ground state, leading to appearance 
of high-momentum and high-energy nucleon distribution in the target. In order to estimate this 
effect in the inclusive cross sections, we consider the phenomenological model [28, 29] where the 
high-momentum (HM) part of the spectral function is determined by excited states with one or 
more nuclei in a continuum. 

In our calculations of the inclusive cross sections only the real part of the optical potential 
is included, because the complex potential produces absorptions of flux. Then, the contribution 
of the Ip- and Is-states to the inclusive cross section (d?a/defdQ.A nrn^jA can be obtained by 
integrating the exclusive cross sections (10) over Q^- The effect of the FSI on the inclusive cross 
section can be evaluated using the ratio 

'^<^""'' = (dr!^) Klfla) ■ P^' 

y^^^fd-^h/RDWIA/ \^^f(^^hJ PWIA 
where (yd^a/defdQ.f) py^jj^ is the result obtained in the PWIA. Then the total inclusive cross section 
can be written as 

d^a ( Sa \ , r ^ -.f d^o- \ ,„^, 



where {d^ a / de jdVi f) h m is the high-momentum component contribution into the inclusive cross 
section [13]. 



III. ANALYSIS OF CCQE INTERACTION AND NEUTRINO ENERGY 

RECONSTRUCTION 

A. Differential cross sections da /d cos 9 and da/dQ"^ 

The charged-current QE events distributions as a function of Q^ or cos 6 were measured by 
K2K [11] and MiniBoone [2, 12] experiments. High statistic data show a disagreement with the 
RFGM prediction. The data samples exhibit significant deficit in the region of low Q^ <0.2 GeV^ 
and small muon scattering angles, which corresponds to forward-going muons. In Refs. [2, 12] it 
was shown that the data/MC disagreement is not due to mis-modeling of the incoming neutrino 
energy spectrum, but due to inaccuracy in the simulation of CCQE interactions. To tune the Fermi 
gas model to the low Q^ region, an additional parameter was introduced, which reduced the phase 
volume of the nucleon Fermi gas at low momentum transfer. In the region of high Q^ the data 
excess is observed, and the values of the axial vector mass Ma, obtained from a fit to the measured 
data, are higher, than the results of previous experiments. 

We calculated the differential cross sections da /d cos 9 and d^a/dQ"^ for neutrino CCQE scat- 
tering off oxygen target in the RDWIA, PWIA and RFGM approaches. 

We note that in the case of (anti)neutrino scattering off free nucleon CCQE the differential 
cross sections [30] da^'^ jdQ'^ at Q^ — > can be written as 

■cos2 0,[F^(O)+Fl(O)] (31) 



da'^''^ G2 



dQ2 27r 
and do not depend on the neutrino energy. The difference 

is proportional to Fa and decreases with neutrino energy. In the range of Ej ~ 0.5 -^ 1 GeV it can 
be used for measuring the axial form factor Fa- 

B. Selection of charged-current QE two-track events 

At the first step, the CC candidate events are selected by requiring that at least one recon- 
structed track must be long and corresponding to a minimum ionizing particle with the momen- 
tum higher, than a few hundred MeVs. The background is originated by neutral-current (NC) 
interactions producing a charged pion. 

In the CC event candidates, the events with one or two reconstructed tracks, with vertex in 
the active target are selected like the CCQE events. No other tracks are allowed to be connected 



10 



with this event vertex. The two-track events are divided into two samples: QE and nQE enriched 
samples. Depending on detector capabilities dE/dx, the information is applied to the second 
track for ir/p separation [31]. Since the QE interaction is a two-particle scattering process, the 
measurement of the muon momentum and angle allows predicting the angle of a recoil proton (the 
second track) assuming the neutrino scattering off to occur with a nucleon at rest. If the measured 
second track agrees with this prediction within A9, it represents likely the CCQE event. Using 
the MC simulation based on the Fermi gas model, the values of A9 are chosen to give a reliable 
separation between the QE and nQE events. 

To study the nuclear-model dependence of this cut, we consider the angle 9pq between the 
direction of outgoing proton and momentum transfer. For neutrino QE scattering off, the nucleon 
is at rest q = Px and cos 6pq =1. For scattering off bound nucleon with momentum p^, it follows 
from Eq.(6a), that 

cos^P, = ^^±^!^. (33) 

The maximum value of 9pq corresponds to scattering off nucleon with a maximum momentum 



Px ~r y Pmax 



COS e™ = ^^ ;^| ~ (34) 



and cos 0™ < cos 9pq < 1. 



In the RFGM the recoil proton energy e^ = \/Pm + ^ — £& + "^ and for \Pmax I — PF we have 



Px=Pf + ^ +'^^yPF + '^'^^ (35) 

where ui = uj — €b- In the RDWIA the energy and momentum of an outgoing nucleon can be written 
(see Eqs.(5), (6a)) as follow: 

Px=Pm + Q (36a) 

Ex = uj + mA — SB- (36b) 

For the scattering off shell nucleon with a maximum momentum p^^ax the energy of recoil nuclei is 



es = \/Pmax +ml^mB+ p^„^/2mij, (37) 

where rriB = mA — 'm + Sm- In the numerical calculations we use |p^(j^|=500 MeV/c and the mean 
missing energy (em)=27.1 MeV for the oxygen target. Using Eqs.(34), (36a), (36b), and (37), we 
have 

cos 9Z = ^i^^ + ^) + iQ'-^')-P'ma^^ (38) 
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where uj = uj — (e^) — Pmax/'^''^*B ^^^ ^B ~ ^^A — in + {e„i)- It follows from (34), that in the 
RDWIA the phase volume in (cos9pg,Q'^) coordinates is larger, than in the Fermi gas model, and 
this difference decreases with momentum transfer. 

C. Reconstruction of neutrino energy 

In the kinematic reconstruction the neutrino energy e^ is formed assuming the target nucleon 
to be at rest inside a nucleus 

_ £/(m - €b) - (el - 2m€b + m^)/2 
[m — eh) — Ef + Kf cos 

This formula ignores the nucleon momentum distribution for the event reconstruction. Using 
Eq.(36a) and the energy balance in the RFGM 



ei + ^/J^^l^m^ - eh = Sf + Ex, (40) 



or 



ei + mA = e^ + Ef + eB- (41) 

in the RDWIA for shell nucleon and 

Ei+EN = ef + ea:. (42) 

for nucleons with energy en in the correlated NN-pair, we obtain the second-order equation for 
the neutrino energy, which takes into account the bound nucleon momentum and the energy dis- 
tributions 

AeI -BEr + C = 0. (43) 

The expressions for coefficients A, B, and C are given in Appendix A for the RFGM and RDWIA. 
The solution of (43) 

Er= (b + Vb^ - 4AC] /2A (44) 

is the reconstructed neutrino energy, which depends on the variables Ip^l, Em, and cost = p ■ 
q/\p- q\- So, the distribution e,.(|p^|,em, cost) corresponds to measured values of (kf, cos 6) and 
at Em , Pm ~^ Eq. (44) has asymptotic form given by Eq. (39) . 
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The n-th moment of Erik f , cos 6 , p„^, Em) distribution versus kf and cosO can be written as 

/Pmax f'£max 

dp / S{p,e)[£r{kf,cos6,p,e)]"-de, (45) 

i-min ^ ^min 

where S{p,£) is the probabihty density function (p.d.f.) for the nucleon momentum and energy, 
the target nucleon momentum and energy distribution being normahzed with respect to the unit 
area. The mean of er{kf,cos6) and its variance o"^(ej.) are defined by 

Erik f, cos 9) = {£r{kf, cos 6)), (46a) 

o-2(e^) = {E^{kf, cosO)) -E^{kf, cos e) (46b) 

In principle, the cut R = a{Er)/Er < S may be imposed (event by event) to select the events with 
well-reconstructed energy. 

The accuracy of reconstructed energy Er{Ei) as a function of Ei can be estimated using the 
moments of Erikf, cos 9) distribution 



(e;^(ei)) = dkf W{kf,cose)[Erikf,cose)]''dcose, (47) 

where W{kf, cos 9) is the p.d.f. of the muon momentum and scattering angle, i.e. 

W{kf,cos9) = \ /" (48) 

ctoti^i) dkfd cos 9 

and 

— — dkf dcos9. (49) 

dkfd cos 9 ■' 

Usually, to select the CC events, kf and cos^ cuts are applied: kf > kcut and cos^ > (cos^)cMi. 
The lower limits of integration in Eq.(47) are {kf)min = kcut, {cos9)min = (cos0)cut and 
[Er{kf, cos 9)]'^ = (e"(fcj,cos0)), if the nucleon Fermi motion effect is taken into account, or 
Er{kf, cos 9) is given by Eq.(39), if this effect is neglected. It's worth to emphasize here, that 
formula (39) can not be used for neutrino energy reconstruction at Ef > (m — ej,) + kf cos 9 or 
Q^ — Qo — 2mej — m'i, because the value of resulting e^ is negative in this region. In terms of 
energy transfer, it corresponds to the range uji < uJ < l02, where uj2 is the solution to the equation 

Ql = 2Ei{Ef - kf cos 9) -ml. (50) 

In the RDWIA, loi is the value of a;, at which 



Ql= \Praax\ + y^^^+^^ " ^^ (51) 
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with Co = uj — {Sm) — Pmaxf^^^^^ ^^^ ^^ ^^^ RFGM, uji is the solution to the equation 

n2 



q: 







Pf + yPp + 2eF0J + oj' 



r,2 



,2 



UJ 



(52) 



where tD = a; — e^ and e\ = p\ + m? . The size of this range Aw = a;2 — wi is proportional to IPmaxI 
{pf) and reduces with increasing cos0. 

The reconstructed neutrino energy e^ = (Sr) is smeared with variance 

a\ei) = {e^,{ei)) - e^.isi) (53) 

and biased with 

A{ei) =ei- Er (54) 

Using this mean energy approach, we estimated the accuracy of the neutrino energy reconstruction 
with and without the nucleon Fermi motion effect in the RDWIA and RFGM approaches. The 
expressions for the moments (e"(A;j,cos0)) and (e"(ej)) are given in Appendix B. 

In the calorimetric reconstruction e^ is formed as a sum of muon energy ej, kinematic proton 
energy Tp and the mean missing energy {em) 

e{kf,cos6) = £f + Tp + (em)- (55) 

The expressions for the moments of £r{kf, cos 9) distribution are given in Appendix B. However, 
the neutrino energy is underestimated in the kinematical and calorimetric reconstructions, when 
the event represents, in fact, the nQE event, but looks like the QE event. 

IV. RESULTS 

The resulting fluxes of neutrino are predicted with the mean energy of --^ 0.7 GeV in the Mini- 
BooNE and T2K experiments, and - 2.5 GeV at the MINOS and MINERvA detectors. We calcu- 
lated the differential inclusive cross sections d^a/dudQ"^ and d? a /dtod cos 6 of CCQE u^ scattering 
off ^^O for these energies using the LEA code, which was adopted for neutrino interaction [13]. In 
Fig.l d'^a/diodQ'^ cross sections, calculated within the RDWIA, PWIA, and RFGM, are shown for 
the neutrino energy ey=Q.7 GeV, and in Fig. 2 for ey=2.b GeV as a function of Q^. Also shown 
in Figures 1 and 2 show the high-momentum component contributions of the nucleon momentum 
distribution in the target. 

At energy ejy=0.7 GeV the RDWIA cross sections, in the maximum, are lower, than the PWIA 
and RFGM results, and this difference decreases with increasing energy transfer. At energy 2.5 
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FIG. 1: (Color online) Inclusive cross section versus the four- momentum transfer Q^ for neutrino scattering 
off ^^O with energy ejx=0.7 GeV and for four values of energy transfer: ^=0.105, 0.159, 0.213 and 0.266 
GeV. The solid line is the RDWIA calculation, whereas the dashed and dash-dotted lines are, respectively, 
the RFGM and PWIA calculations. The dotted line is the high-momentum component contribution to the 
inclusive cross section. 



GeV, in the range of lj > 0.5 GeV, around the peak, the RFGM results are lower, than the 
RDWIA ones and fall down rapidly as Q^ decreases. This trend is characteristic of the nucleon 
momentum distribution and the Pauli blocking effect as calculated in the Fermi gas model. On 
the contrary, the Q^-dependence of the RDWIA and PWIA cross sections at low Q^ is softer due 
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FIG. 2: (Color online) Same as Fig.l, but for the neutrino energy £,^=2.5 GeV and for four values of energy 
transfer: uj=0.279, 0.507, 0.735 and 0.962 GeV. 



to the HM-component contribution, which becomes dominant at Q^ < 0.1 (GeV/c)^. 

Generally, theoretical uncertainties of the correlated NN-pairs contribution to the inclusive cross 
sections are higher as compared to the shell-nucleons contribution. The electron-nucleus scattering 
data [14, 32, 33] show that more complicated configurations, than a simple hard interaction between 
two nucleons, are involved in this case. Moreover, the off-shell ambiguities will be important for 
the high-momentum component, and one might expect the details of the off-shell extrapolation to 



16 




> 

(D 
O 



o 

X) 



b 

X) 



o 



16 

14 

12 

10 

8 

6 

4 

2 



0,2 0.4 0.6 

Q' [GeV/c]' 



- 












^ 


r" ^ ■ 


E.= 


=1.2 GeV 




- 


/ 
/ 




\ • 






: 


"1 




\ \ ^. 






— 


-1 ■ 

1 

-.' 

-r 

f 




\ ^ 
\ s 




^ 


^ 


-/.[ 













> 

o 



o 

XI 



b 

CO 

o 



> 

cu 
O 



O 

X) 



^b 

XJ 



O 




0,2 0.4 0.6 0. 
Q' [GeV/cf 



0.2 0.4 0.6 0. 
Q^ [GeV/c]" 




0.2 0.4 0.6 0. 
Q' [GeV/cf 



FIG. 3: (Color online) Inclusive cross section versus the four-momentum transfer Q^ for neutrino scattering 
off ^^O and for four values of incoming neutrino energy: e^=0.5, 0.7, 1.2 and 2.5 GeV. The solid line is 
the RDWIA calculation, whereas the dashed and dash-dotted lines are, respectively, the RFGM and PWIA 
calculations. The dotted line is the high- momentum component contribution to the inclusive cross section. 



become critical [34]. 

The inclusive cross sections for energies £,^=0.5, 0.7, 1.2, and 2.5 GeV are presented in Fig. 3, 
which shows cfa/Q'^ as a function of Q^. Here the results, obtained in the RDWIA, are compared 
with cross sections calculated in the PWIA and RFGM. The contributions of the NN-correlations 
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FIG. 4: (Color online) Inclusive cross section versus the niuon scattering angle for four values of incoming 
neutrino energy: e,j=0.5, 0.7, 1.2 and 2.5 GeV. The solid line is the RDWIA calculation, whereas the dashed 
and dash-dotted lines are, respectively, the RFGM and PWIA calculations. 

are shown as well. The cross sections, calculated in the Fermi gas model, are higher, than those 
obtained within the RDWIA, and this difference increases with decreasing Q^. At Q^=0.1 (GeV/c)^ 
this discrepancy equals 54% for £,^=0.5 GeV and 43% for £,^=2.5 GeV. In the region around the 
maximum (5^=0.2 (GeV/c)^ the difference is about ~18% for 8^=0.5 GeV and --^11% for £,^=2.5 
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GeV. At (5^=0.05 (GeV/c)^ the contribution of the HM-component increases with energy from ~ 
15% up to 23% in the energy range of 0.5 -i- 2.5 GeV. 

Figure 4 shows the inclusive cross sections da /d cos 6 calculated in the RDWIA, PWIA, and 
RFGM approaches for energies £,^=0.5, 0.7, 1.2, and 2.5 GeV. They are displayed as a function of 
cos^. It is clear that in the region 0.8< cos^ <1 the values of the RFGM cross sections are higher, 
than those obtained within the RDWIA, and this difference decreases with neutrino energy. For 
energy e,^=0.5 GeV (e^=2.5 GeV) this discrepancy is about 25 times (~11%) at cos 0=0.95 and 
~89% (~2%) at cos 0=0.8. We note that measured Q^ and cos ^-distributions of the CCQE events 
[2, 11, 12] show similar features as compared to the RFGM prediction. 

Figures 5 and 6 show cos 0™ as a function of Q^ calculated within the RDWIA and Fermi gas 
model kinematics for energies £1^=0.5 and 2.5 GeV. The outgoing proton carries the kinematic 
energy, that is approximately to. So far as to is low, the problem consists in identifying the events 
with very soft recoil proton; for high uj this proton has high energy and may interact in the detector, 
making particle identification and track reconstruction more challenging. In these figures we show 
the contours of the phase volume in the {cos9pq,Q'^) coordinates for 0.25< uj <1 GeV. Apparently, 
in the RDWIA kinematics this volume is larger, than in the RFGM. On the other hand, the 
difference decreases with uj and neutrino energy. Thus, systematic errors for the efficiency and 
purity of the two-track events selection are nuclear-model dependent. 

We have studied the accuracy of the neutrino energy reconstruction with neglecting the sys- 
tematics related to the event selection and resolution, i.e. with no detector effects or background. 
The study was performed with the values of cuts {kf)cut=0.2 GeV/c and (cos0)cut=O. 

In Fig. 7 the uncertainties of the energy reconstruction using Eq.(39) within the RDWIA and 
RFGM approaches are presented as functions of neutrino energy. The top panel shows the bias 
A = (£j— £r)/£j, the middle panel shows the variance a/ei (the energy resolution), and the efficiency 
of the one-track events detection is displayed in the bottom panel. It is clear that in the case of 
the Fermi gas model Eq.(39) systematically underestimates the neutrino energy and A decreases 
as the energy increases from —4.7% for £,^=0.3 GeV up to —0.7% for £,^=2.5 GeV. The variance 
a/si (efficiency) increases with energy from --^5.4% (~71%) up to -^12% (~99%) over the range of 
energy from 0.3 to 2.5 GeV. 

In the RDWIA approach A r^ -5.2% for £^=0.3 GeV and A ?«3.9% for £i.=1.6 GeV. At energies 
£jy >1.6 GeV and at fixed values of cuts the denominator in Eq.(39) can be negative in the detection 
volume of the {kf, cos9) phase space, which is determined by Eqs.(50) and (51). We note that 
the size of this volume decreases as {cos 9)cut increases. The energy resolution is about 8.3% for 
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FIG. 5: (Color online) Contours of the phase volume in the (cos^pg, Q^) coordinates for neutrino scattering 
off ^^O with energy £1^=0.7 GeV and for four values of energy transfer: Ci;=0.288, 0.320, 0.374 and 0.427 
GeV. The solid line is the RDWIA calculation, whereas the dashed line is the RFGM calculation. 



eu=0.3 GeV and ~12.7% for ey=l.Q GeV, and the maximum of 15.4% is located around ei^=0.8 
GeV. The efficiency rapidly increases with energy from 28% for £,^=0.3 GeV up to 96% for £u=l.Q 
GeV. So, the values of the bias and energy resolution, obtained within the RDWIA, are higher, 
than those obtained in the Fermi gas model. 

The accuracy of the mean energy method, which takes into account the nucleon momentum 
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FIG. 6: (Color online) Same as Fig. 5, but for neutrino energy £,^=2.5 GeV and for four values of energy 
transfer: w=0.279, 0.507, 0.735 and 0.962 GeV. 



distribution in the target, is shown in Fig. 8. We assume that the maximum neutrino energy in 
the neutrino beam (Eq.(B8)) is Efnax=W GeV. The top and bottom panels show the biases and 
energy resolutions calculated within the RDWIA and RFGM. In the case of the RFGM the mean 
energy method systematically overestimates the neutrino energy; A=4.3% for 8^=0-3 GeV, A «6% 
for £,^=0.8 GeV and decreases down to 3% for ei,=2.5 GeV. The energy resolution increases with 
energy from 4.6% up to 11% in this energy range. In the RDWIA approach A ~ —4%, {cr/ei) ~14% 
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the efSciency (bottom panel) of the one-track events detection with kf >0.2 (GeV/c) and cos 6* >0 as 
functions of neutrino energy. The neutrino energy reconstruction was formed assuming the target nucleon 
to be at rest inside the nucleus. The vertical bars show (T[{ei — er)/£i]- As displayed in the key, the biases, 
variances and efficiencies were calculated in the RDWIA and RFGM. 

for e,y=0.3 GeV, and A rs4.5%, {cr/si) f«10.5% for £,^=2.5 GeV. It should be noted here that bias 
may depend on the value of Emax- 

The effect of the nucleon momentum distribution in the target is shown in Fig. 9. The biases, 
calculated within the RDWIA (top panel) and RFGM (bottom panel) using Eq.(39) ( Aj^ ) and 
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the mean energy method ( A^e )) are presented as functions of neutrino energy. In the RDWIA 
approach the Afr and A^e show similar behavior with neutrino energy, and the nucleon Fermi 
motion effect leads to increasing the bias by about 1.2%. In the Fermi gas model with this effect e^ 
is overestimated, and A/^,(A„e) = -4.7%(4.3%) for energy 0.3 GeV and Afr{Ame) = -0.7%(3.4%) 
for e„=2.5 GeV. 
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Apparently, the accuracy of the kinematic reconstruction of neutrino energy for one-track 
events depends on the nuclear models of QE neutrino CC interaction with nuclei and on the 
neutrino energy reconstruction methods. In the K2K and MiniBooNE experiments Eq.(39) is 
applied for the energy reconstruction. The bias Apc and energy resolution 6fg = {'^/^i)FG 
were calculated using MC simulation based on the Fermi gas model. We can estimate the 
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FIG. 10: (Color online) Bias (top panel), variance (middle panel) of the reconstructed neutrino energy, and 
the efficiency (bottom panel) of the two-track events detection with kf >0.2 (GeV/c) and cos9 >0. The 
vertical bars are the same as in Fig. 7. The bias, variance and efficiencies were calculated within the RDWIA. 

systematic uncertainties of this approach by comparing Apc and 6fg with Ar and 5r eval- 
uated in the RWDIA approach using the mean energy method. It is clear that uncertainties 
depend on neutrino energy: ApGi^R) ~ -4.7%(-4%) and dpciSR) ^5A%{13.7%) for e^=0.3 
GeV; Afg{Ar) « -2.3%(4.1%) and 6fg{Sr) f=^10.6%(16.3%) for e^=0.8 GeV; Afg(.Ar) « 
— 0.7%(4.5%) and 6fg{^r) ~12%(11.5%) for £,^=2.5 GeV. So, the bias uncertainty increases with 
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energy from (A/j — A^^) RiO.7% for e^=0.3 GeV up to 5.2% for e,^=2.5 GeV, and the energy 
resolution uncertainty decreases with increasing energy from 6^ — 6fg ~8.3% down to 0.5% in this 
energy range. We note that these estimations may depend on the values of {kf)cut, {cos 6)cut and 

In Fig. 10 the accuracy of the energy reconstruction for the two-track events, calculated using 
Eq.(55) within the RDWIA, is shown as a function of neutrino energy. The top panel shows the 
bias, the middle panel shows the variance, and the bottom panel shows the efficiency of the two- 
track events detection with cuts kf >0.2 GeV/c and cos^ >0 for muon tracks and without any 
cuts for the proton tracks. At energy Si, >0.3 GeV the bias is A ~ —0.1% and does not depend on 
the neutrino energy. The energy resolution decreases as the energy increases from 3.4% for £,^=0.3 
GeV up to 0.5% for ejy=2.5 GeV, and the efficiency rapidly increases with energy from ~16% up to 
~44% in this energy range. These estimations show that the accuracy of the calorimetric method 
can be higher, than the kinematic one and does not depend on the model of CC neutrino QE 
interaction with nuclei and on the nucleon momentum distribution in the target. The challenge is 
identifying the proton track and reconstructing its kinetic energy with reliable accuracy at the low 
threshold energy for proton detection. 

V. CONCLUSIONS 

In this paper, we study the quasi-elastic neutrino charged-current scattering on the oxygen 
target in various approximations (PWIA, RDWIA, RFGM) making particular emphasis on the 
nuclear-model dependence of the results. In the RDWIA, the LEA program, adapted to neutrino 
interactions, was used for calculating the differential cross sections with the effect of NN-correlations 
in the target ground state. 

The inclusive cPa/dQ'^ and da /d cos 6 cross sections, calculated within the RDWIA, and the 
measured Q^, cos ^-distributions of CCQE events exhibit similar feature as compared to the Fermi 
gas model. The magnitude of inclusive cross sections SajdQ'^ and da /d cos 9 is lower in the 
RDWIA calculations, than that of the Fermi gas model, and in the region around the maximum 
Q2=o.2 (GeV/c)2 the difference is about 18% for e^=0.5 GeV and 11% for e^=2.5 GeV. The 
contribution of the HM-component at (5^=0.05 (GeV/c)^ increases with neutrino energy from 15% 
up to 23% in this energy range. Note that the measured Q^ and cos ^-distributions of CCQE events 
are also lower, than the RFGM prediction at low Q^. 

We have shown that the efficiency and purity of the CCQE two-track events selection are nuclear- 
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model dependent, and the difference decreases with increasing energy transfer and neutrino energy. 

We have studied the nuclear-model dependence of the energy reconstruction accuracy, neglect- 
ing the systematics related to event selection and resolution. We found that the accuracy of the 
kinematic reconstruction for one-track events depends on the nuclear model of CCQE neutrino 
interaction and on the neutrino energy reconstruction method. The uncertainties in the recon- 
structed energy bias increase in the energy range of 0.3 -^ 2.5 GeV from ~0.7% up to 5.4%, and 
the energy resolution ambiguities decrease from 8.3% down to 0.5% with increasing energy. In the 
case of two-track events the accuracy may be higher and does not depend on the nuclear models 
of CCQE neutrino-nucleus interaction. 

We conclude that the use of RDWIA in the Monte Carlo simulation of neutrino detector and 
the data analysis would allow one to reduce the systematic uncertainty in neutrino oscillation 
parameters. 
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APPENDIX A: EQUATION FOR NEUTRINO ENERGY 

In Eq.(43) 

Ae^ -Bsr + C = 0. 

the coefficients A, B, and C are defined as follows: 

A = a^ - pI^z"^ , (Ala) 

B = ab-2pl,z^kf cose, (Alb) 

C = b^/A-plz^k}, (Ale) 

where 

a = Eef — kf cos 6, (A2a) 

b = elf-{\pj^ + m^)-k}, (A2b) 

z = cosT = p^-q/\p^-q\. (A2c) 
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In the Fermi gas model 



Eef = Sf - v^p2^ + m2 + €b, (A3) 

in the RDWIA model, for shell-nucleons 

Eef = ef-mA + [{mA -m + e^f + P^]^^^ (A4) 

and for nucleons in the correlated NN-pair 

APPENDIX B: MOMENTS OF RECONSTRUCTED NEUTRINO ENERGY 



In the Fermi gas model with the p.d.f. 



3 



S{Pm: ^m) = T^^ni - Cfe) (Bl) 



Eq.(45) takes the form 



O rPmax rZmax 

(e"(/cf, cos 6*)) = TT I p^dp I [e.r{kf, cos 9, p,z)Y^dz, (B2) 

47rni y„ . J , - 

where z = cost and e^ is given by Eqs.(44), (Al), (A2), and (A3). In the RDWIA the p.d.f. can 
be written as follows: 

S[v, s) = X] ^"'S'a(P)'^[^ ~ (em)^] + vnnSnn{p, e) (B3) 

a 

and we have 

{el' {kf, cos 6)) = ^VaiSr^ai^f, cos d)) + VNN{£r,NN{kf, COS 6)) , (B4a) 



a 

Pmax 



{e'^. ^^{kf, COS &))= p^dp 5„(p)[e,.,o(/c/,cose,p,z)]"dz, (B4b) 



mtn ^ '^rrnn 



/Pmax r^max P^f-nax 

p^dp I del SNN{p,£)['^r,NN{kf,cos9,p,e,z)]'^dz, (B4c) 

where 5*0, and Snn are, respectively, the p.d.f. for the momentum and energy of nucleons on the 
shell a and in the correlated NN-pairs, e"^ and e"^ ^^ are given by Eqs.(44), (Al), (A2), (A4), 
(A5) and the sum is taken over occupied shells. The coefficients Va and vnn are 

Va,NN = {-jr~i a) ' (B5) 

0'a,NN ^ dkfd COS t> / a,NN 
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where 

/,2 
(37""} o) dkfd COS 9 (B6) 

KdkfdcosOJ a,NN ■' 

The integral is calculated with {kf)min = {kf)cut and {cos9)„Lin = {cos9)cut- Using Eqs.(47), (B4), 
we have 

(e"(ei)) = X] ^"(^r (£i))a + w'ArAr(e"(ei))iVAr, (B7) 

where Wa,NN = (^a,NN /{Yla ^a + o'nn)- In Eqs.(B2) and (B4) the limits of integration over z are: 
-2mm=-l and Zmax=min{l, z/}. The value of z; is obtained from the requirement er{kf, cos 9) < 
Emax, where Emax is the maximum neutrino energy in the neutrino beam. We note that this 
constrain on Zmax leads to increasing bias A in the reconstructed energy. Then from Eqs.(36a), 
(40), (41), and (42) with Cj = Emax it follows, that 

where ujmax = Emax — £/ and q^ = E^ax + k'j — 2Emaxkf cos 9. In the Fermi gas model e = 



\/Pm + "i^ — Cfe) in the RDWIA for scattering off shell-nucleons e = niA — \/Pm + "^B' ^^'^ ^ = ^N 
for scattering off nucleons in the correlated NN-pair. 

For the two-track events the moments of the er{kf,cos9) distribution can be written as 

{e'^{ei)) = Y,^e.{erie^))a, (B9) 

a 

where 



'max 



(e;!(ei)) = I dkf I dcos9 / dcp [ej + Tp + {em)rWa{kf, cos 9,(1), pm)dpm, (BlOa) 

JO ^ Pmin 



^ a 



d^a 



(BlOb) 



dkfd cos 9 dcjidpn 
dkf I dcos9 I d<p I dk fd cos 9d<pdp. -^^-' ^^'°') 



^1 UK,f I U UOS U 

/O Jp^ 



mm 



Wn = cr. 



7/Y.< (BlOd) 



and d^ a /dkfd cos 9 d(j)dpm is the QE neutrino CC scattering exclusive cross section (10) in terms of 
a missed momentum pm- 
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